Identification of Non-unitary triplet pairing in a heavy Fermion superconductor UPt3 
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A NMR experiment recently done by Tou et al. on a heavy Fermion superconductor UPt3 is 
interpreted in terms of a non-unitary spin-triplet pairing state which we have been advocating. 
The proposed state successfully explains various aspects of the seemingly complicated Knight shift 
behaviors probed for major orientations, including a remarkable d-vector rotation under weak fields. 
This entitles UPt3 as the first example that a charged many body system forms a spin-triplet odd- 
parity pairing at low temperatures and demonstrates unambiguously that the putative spin-orbit 
coupling in UPt3 is weak. 
PACS numbers: 74.70.Tx, 71.28.+d, 74.60.-w 



Except for the famous example of neutral superfluid 
Fermion system 3 He, no charged Fermion many body 
system has been known to exhibit spin triplet pairing 
in nature so far; that is, a triplet superconductor is not 
firmly identified yet. A heavy Fermion superconductor 
UPt3 [|) has been regarded as a prime candidate for a 
triplet superconductor because the observed double tran- 
sition (T c i ^0.58K and T c2 ~0.53K) and multiple phases 
(A, B, and C phases) in H vs T apparently point to some 
unconventional pairing state where a Cooper pair must 
have internal degrees of freedom responsible for the mul- 
tiple phases. 

A central controversy in identifying the pairing symme- 
try in UPt3 lies in whether the internal degrees of freedom 
as the origin nearly degenerate superconducting transi- 
tion temperatures come either from the orbital part [^]||] 
or from the spin part 0] of a Cooper pair. There are the 
three major scenarios for explaining the multiple phase 
diagram, apart from the so-called accidental degeneracy 
scenario || whose proposed state is not qualified in detail 
as yet, being unable to compare it with the above three on 
an equal footing. The orbital scenarios identify the pair- 
ing symmetry as a two-dimensional representation whose 
multi-dimensionality of the orbital part is responsible for 
the degeneracy, namely the singlet pairing E l£( proposed 
by Joynt ||, the triplet pairing E 2u by Sauls The 
latter assumes that the degeneracy is broken by the anti- 
ferromagnetism and spin-orbit coupling (SOC) is strong. 
Thus in F,2 U the spin of a pair is firmly locked to the crys- 
tal lattice. In the spin scenario j|, on the other hand, 
the degeneracy comes from the spin part of the pair in 
which weak SOC is assumed by considering that SOC 
in a one-particle level is already taken account in form- 
ing heavily renormalized quasi-particle characteristic of 
heavy Fermion materials. 

It was not an easy task to distinguish the above 
three scenarios. Since usual thermodynamic or trans- 
port measurements extract power low indices in their In- 
dependence, which reflect the energy gap topology. The 
accumulated experiments all point to the gap vanish- 



ing at both point and line This is consistent with 
the above scenarios which yield more or less similar gap 
topology. An obvious and decisive experiment is to mea- 
sure the Pauli spin susceptibilities Xi f° r the principal ori- 
entations (i — a,b, c-axes) in the hexagonal crystal since 
in the singlet Ei 9 , Xi decreases for all directions, and in 
the triplet E2 U , Xa and Xh (Xc) remain unchanged (de- 
creases) below T c , irrespective of the three phases (A, B, 
and C), that is, the three phases give rise to the identical 
behavior in their scenarios. In the spin scenario the Xi 
behavior is much more subtle and complicated as shown 
shortly. 

Owing to technological developments in making high 
quality single crystals by Onuki group and improvements 
of NMR instrumentation by Asayama and Kitaoka group, 
their joint work ^| has culminated in success of measur- 
ing the Knight shift (KS) by probing 195 Pt nuclear mag- 
netic resonance using single crystals (two crystals: |J3, 
jJ4). They extend and refine the previous experiment Q| 
to much more wider fields and yet unexplored field ori- 
entations. The new experiment covers virtually all the 
phases in H vs T plane for major H orientations, en- 
abling us to completely determine the Cooper pair spin 
structure. 

The purpose of this paper is to demonstrate that a sin- 
gle and simple Ginzburg-Landau (GL) free enegy based 
on a few basic assumptions can explain all aspects of 
the seemingly complicated results of their experiment j^] 
summarized in Table 1, leading us to identify the pairing 
symmetry in UPt3. Before going into detailed analysis of 
the Knight shift changes for major orientations, we first 
reexamine our fundamental assumptions in the light of 
the new NMR experiment. 

The first and most central assumption, that is, the ef- 
fective SOC felt by Cooper pairs is relatively weak. This 
implys that under symmetry operations in D^h x 5*0(3) x 
T x U(l) (Deh hexagonal point group, SO(3) spin rota- 
tion, T time reversal and f (1) gauge symmetry) the spin 
and orbital parts of the spin-triplet order parameter de- 
scribed by A Q ,p(k) — i(cr ■ d(k)<72)a,;8 in terms of the d 
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vector behave independently. The strong SOC assump- 
tion adopted by the orbital scenarios, under the symme- 
try operations the spin and orbital part transform simul- 
taneously and the Cooper pair spin is strongly locked to 
the crystal lattice and is not rotatable. In this strong 
SOC case one must resort to the remaining possibility of 
the orbital part, which is the case in the F,2 U scenario. 
Tou's result @ that the field as low as ff(|| c) ~2kG 
makes the Cooper pair spin rotate clearly shows that 
weak SOC is indeed the case for UPt3. (Note that the 
BC boundary is far off at ~12kG.) This is one of the 
most important significance in their experiment because 
almost all theories for heavy Fermion superconductors 
take strong SOC for granted since Anderson's proposal 

The second assumption of spin-triplet degeneracy of 
our pairing function is consistent with Tou' experiment 
H since through their Knight shift change one can esti- 
mate the Pauli spin susceptibilities Xi f° r various direc- 
tions coming from the quasi-particles at the Fermi sur- 
face. (Don't confuse those with the total bulk suscep- 
tibilities which are of course very anisotropic as is well 
known) According to their data 0, Xi both parallel and 
perpendicular to the c-axis are nearly isotropic, support- 
ing the triple degeneracy in spin space and leading to 
nearly SO(3) symmetry. 

The third assumption that the antiferromagnetic (AF) 
fluctuations characterized by the triple-Q vectors with 
Qi = (^,0,0) and its equivalent positions Q2 and Q3 in 
reciprocal space (a* = b* = -^=, c* = 2 ^ L ) in the hexag- 
onal plane pd| ] are responsible for breaking the triple de- 
generacy. This is not inconsistent with their experiment 
because the static AF order is not present M . The triple 
Q-AF fluctuation structure probed by neutron experi- 
ments H |ll| is symmetry-broken from the outset by, for 
example, the incommensurate lattice modulation whose 
presence or absence exhibit a good correlation with the 
double transition |l2| , leading to a different weight to the 
triple Q fluctuations. Therefore it is possible to break 
the hexagonal symmetry by this longitudinal AF fluctu- 
ations. Without loosing generality we take Qi \\ a* || b 
whose fluctuation is differed from those at Q2 and Q 3 , 
implying that the symmetry of the hexagonal plane is 
lowered to two-fold symmetry, which is indeed observed 
by Tou's experiment |(| as will be seen. They single 
out the two-fold unique axis in the nominally hexagonal 
plane. That is remarkable. 

Let us now start out by writing down the phenomeno- 
logical GL free energy, derived based on the above three 
assumptions: 

F = ao(T-T c0 )|d| 3 +i J 9 1 |d|*+i/3a|d 2 | a -7|b.d| a -A|c-d| a 

where Qi || b. The last term A|c-d| 2 which is somewhat 
ad hoc expresses a weak anisotropy of the order parame- 



ter in spin space, reflecting the fact that the Knight shift 
changes below T c for H=2kG parallel to the c-axis [||. 
This free energy leads to 

f= « (r - T^)|^] 2 + |/3i|d| 4 + i^ld 2 ! 2 

j—a,b,c 

withT c b = T c0 + J- > T c c = T cQ + ± > T- for 7 > A > 0. 
The three transition temperatures T", Tj? and T£ which 
were originally degenerate are split into the two groups: 
Tj? and {T° and T£} by the symmetry breaking field as 
mentioned. The latter two are further assumed to be 
slightly different T® < T c c due to the A-term because of 
the small uniaxial anisotropy of the system. 

At T = T c h , which is identified as the upper critical 
temperature T c \ = 0.58/^ the A phase characterized by 
d a = 0, rfh ^ and d c = appears first. Then at a 
lower temperature, T — Tj? which is identified as T C 2 — 
0.53 A" the second order transition from the A phase to 
the B phase characterized by a non-unitary triplet state: 
d a = 0, db 7^ and d c ^ takes place where the above 
free energy leads to the phase difference between db and 
d c by 7r/2 when > 0. It can be proved within the 
above free energy that third transition at T — T c a , which 
is designed to situate a few mK below T — T c c never 
realized at zero field. 

Let us now discuss the phase diagrams for external 
field H applied parallel to and perpendicular to the o 
axis. We must take into account the vortex structure 
based on the above free enegy functional by adding the 
terms describing the spatially varied order parameter: 
F = Fg rai i + Fbuik with 

F grad = i K i(\ D * d i\ 2 + lA^f) + Ki\D z dj\ 2 } 

j=a,b,c 

F b uik= £ a {T -Tl)\dtf + \f3 1 \d\ i 

j—a,b,c 

+i/3 2 |d 2 | 2 + iA X p|H.d| 2 . 

The gradient term F grac i describes the spatial variation of 
the order parameter under external field with the gauge 
invariant derivative Dj — —ihdj — —Aj (A is the vector 
potential). Here there exist two kinds (K\ and K ] 2 (j = 
a, b and c) of the gradient term associated with hexag- 
onal symmetry where K\ {K%) can differ from Kf and 
Kl (K2 and iff) because of the symmetry breaking AF 
fluctuation, that is, K\ =K X - CM 2 , K\ = K 2 - C'M 2 , 
Kf = A1+CM 2 , K% = A 2 +C'M 2 , Kl = K[+(M 2 , K% = 
K' 2 + C'M 2 . Here M 2 is the amplitude of the AF fluctu- 
ation polarized along the 6-axis. It should be noted that 
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there is no so-called gradient mixing term, which washes 
out the desired tetra-critical point, in the present spin 
scenario. This is quite different from that in the orbital 
scenarios ^|,[| where the gradient mixing is inevitable. 

The resulting phases in the H vs T shown in Fig.l 
are fully characterized by specifying the d vector in each 
phase, that is, its direction, the number of the compo- 
nents, and the relative phase. The three vectors in the 
spin-triplet state are now denoted by the real vectors a, b 
and c referring to the Cartesian coordinate in the hexag- 
onal crystal. The high T and low (high) field A(C) phase 
is described by a single component d vector, while the 
B phase in the low T and low H is a non-unitary state 
characterized by two-component d vector whose relative 
phase is 5. Because of the Axp term in the above free 
energy, which tends to align the d- vector perpendicularly 
to H, the B phase for H || c is further subdivided into 
b+ic in low H and b+ia in high H . As H(\\ c) increases 
the d vector rotates, corresponding to the KS behavior 
at H -2kG in Table 1. 

We analyze the KS experiment by Tou et al. HQ, 
whose results are summarized in Table I: In the case 
of H || c the KS does not (do) change below T c for 
H < 2kG (H > 3kG) (see their Fig.4 §). This charac- 
teristic field H rot differs from that of the BC transition at 
12kG, through which KS remains unchanged. According 
to the above theory within the B phase region b+ic in 
low H changes into b+ia in high H because the d vector 
tends to orient perpendicularly to the field direction, and 
the BC transition is from b+ia to a. Thus KS change 
does occur at H rot , never at Hbc- This is precisely the 
case seen experimentally. The implication of the d vec- 
tor rotation is obvious that SOC to lock the spin to the 
lattice is not strong, but weak enough to reoriente under 
a field as low as 2kG. 

When H || a, since both b+ic (B phase) to c (C phase) 
are perpendicular to H, KS should not change below T c 
for any field strengths, even crossing the BC boundary 
situated at ^7kG. This coincides with the observations 
done for #=11, 4.7, 2.1 and 1.8KG (see their Fig.2 §]). 

As for H || b, KS should and should not change below 
T c for the B phase and C phase respectively. The exper- 
iments under H =1.9, 2.4, 3.4, 4.7 and 5.3kG (B phase) 
and H =8.0 and 9.7kG (C phase) confirm this prediction 
(see their Fig. 3 This means that the SOC in the 

basal plane is strong enough not to depin the d vector 
parallel to the the &-axis by the field up to Hbc ^7kG. In 
view of the difference in the characteristic energy scales 
for H || b (T c b -T c c -50mK) and H || c (T r c -T L a <50mK), 
this is understandable. 

As pointed out by Tou et al. S, their data allow to 
draw the information on the A phase, which was quite 
scarce. By carefuly comparing the two lowest field data 
sets for H \\ c (their Fig.4) and H \\ b (their Fig. 3) one 
can clearly see that the decrease of KS does begin to start 



at T cl (T c2 ) for H \\ b (H || c), indicating that in the A 
phase the d vector points parallel to the 6-axis. This is 
precisely consistent with our theory. 

Importantly enough, the basal plane shows the two- 
fold symmetry The 45° direction data from the a-axis 
exhibits an intermediate decrease in the KS change. The- 
oretically we expect that the decrease in x corresponds 
to cos 2 (H • d) = \. 

The absolute values of the KS change (~0.06% (||b) 
vs ~0.07% (||c)) which are proportional to the Pauli 
spin susceptibilities Xi coming from the quasi-particles 
near the Fermi surface are roughly isotropic when taking 
into account the anisotropic hyperfine coupling constants 
(~ 20% difference), supporting the spin space degeneracy 
scenario (see their detailed analysis ||). The magnitude 
of the estimated Xi is extremely small; only a few % of 
the total KS, but Tou's experimental accuracy can sen- 
sitively detect it, negating a persistent opinion that the 
KS measurement is meaningless because the susceptibil- 
ity is dominated by the T-independent van Vleck term. 
It is indeed dominated by the van Vleck term, but they 
are undoubtedly probing the relevant quasi-particle con- 
tribution or the Pauli susceptibility part. 

An immediate consequence of this identification al- 
lows us to determine the direction of the spontaneously 
induced moment M s associated with the non-unitarity 
in the B phase: M s a id x d* || a at low fields where 
d = b+ic. This spontaneous moment M s was observed 
and its magnitude was estimated as ~0.5G by /iSR exper- 
iment Now the direction of M s should be probed. 
Since macroscopically this moment is cancelled out by 
the screening surface current and thus the magnetic in- 
duction B = in the bulk, the unscreened M s survives 
only near the surface whose depth is an order of the pen- 
etration length. If the supercurrent whose velocity is v s 
(say, lem/sec) flows along the c-axis, the electric field 
E = v s x M s is spontaneously induced transversely, i.e., 
across the 6-axis, whose order of magnitude is estimated 
as ~ 10" 7 volt/m. This non-trivial prediction should be 
checked experimentally. 

Having analyzed Tou's experiment we discuss 

some of the critiques to our theory in the light of the 
new findings. The crossing phenomenon of the two up- 
per critical field curves for H c2 _L c and || c has been 
taken as strong evidence that the Pauli limiting is absent 
(present) for the former (latter), supporting the E2 U sce- 
nario (m). Tou's experiment |(|[7j shows that this is not 
the case because both %'s remain unchanged below T c . 

Although we have determined the spin structure of the 
Cooper pair fairly completely, the precise gap topology 
remains to be seen. Recent several direction-resolved 
transport |l5j and thermodynamic Jl6| measurements 
for each phase should be reexamined in the light of the 
present theory. We emphasize that the present scenario 
gives rise to different gap structures for the B, and A 
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and C phases because the former (latter two) state is 
non-unitary (unitary), which have distinctively different 
low enegy excitation spectral structures. 

In conclusion, we have identified the spin structures of 
the Cooper pair in each phase of UPt3, which are fully 
consistent with the NMR experiment by Tou et al. Q . It 
is rather self-evident now that since the KS data show at 
least the four different behaviors for the principal three 
directions, it is impossible to understand in terms of the 
two-component pairing state such as the orbital scenar- 
ios. The main attractive force for the spin-triplet pairing 
comes from the paramagnetic fluctuations observed be- 
fore 

[0- The triple de generacy is weakly broken by the 
AF fluctuation, resulting in one of the complex d vec- 
tor component locked to the 6-axis. The fact that the 
remaining component of the d vector is rotatable under 
weak fields verifies experimentally that the effective spin- 
orbit coupling felt by Cooper pairs is weak. 

The authors thank H. Tou, Y. Kitaoka, Y. Onuki, T. 
Sakakibara and A. Sawada for useful experimental infor- 
mation. In particular, we are much indebted to H. Tou 
who has performed successfully this difficult NMR exper- 
iment and to his day-to-day information. 
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TABLE I. Summary of the Knight shift experiment for the 
principal directions a, b, and c-axes. x (o) denotes the Knight 
shift unchanges (changes) below T c . 
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FIG. 1. Schematic phase diagrams in H vs T for H \\ a, b 
and H || c where the direction and the components of the d 
vector for each phase are indicated. (Hbc the second phase 
transition between the B and C phases, H rot the rotation field 
of the d vector) 
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